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The maximum efficiency of autonomous engines with finite chemical potential difference is inves¬ 
tigated. We show that without a particular type of singularity autonomous engines cannot attain 
the Carnot efficiency. This singularity is realized in two ways: single particle transports and the 
thermodynamic limit. We demonstrate that both of two ways really lead to the Carnot efficiency 
in concrete setups. Our results clearly illustrate that the singularity plays a crucial role for the 
maximum efficiency of autonomous engines. 
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Introduction. —The maximum efficiency of heat en¬ 
gines has been one of the central issues in thermody¬ 
namics. Carnot showed that the efficiency of an engine 
attached to heat baths with temperature Th and 
{Th > Tif) is bounded by 1 — Tl/Th [i|. The upper 
bound is attained when the external control on the en¬ 
gine is quasistatic. For the case with two particle baths 
under isothermal condition, the efficiency is bounded by 
1. These maximum efficiencies are called the Carnot ef¬ 
ficiency (CE). Nonequilibrium thermodynamics has re¬ 
cently been applied to small fluctuating systems with ex¬ 
ternal control, where the maximum efficiency analogous 
to macroscopic thermodynamics has been established . 


Since most engines from electric power plants to molec¬ 
ular motors are autonomous, thermodynamics for au¬ 
tonomous engines is an important issue. Here, the word 
autonomous stands for systems with no time-dependent 
control parameter in non-equilibrium steady states Q. 
Since all variables of the engine inevitably fluctuate, the 
attainability of the CE with autonomous engines is a 
non-trivial problem. In the linear response regime (i.e., 
Th — Tl — 0), it is well known that the tight-coupling 
condition is necessary for autonomous engines to attain 
the CE [J. On the other hand, for the case of finite 
difference of temperatures or chemical potentials, most 
of studies have paid attention to specific and elaborated 
models Famous models of autonomous engines 

are Feynman’s ratchet Q and the Biittiker-Landauer sys¬ 
tem [y, Q , which convert heat flux into work. Although 
they seemingly attain the CE [ 1 , § , it has been estab¬ 
lished that these models actually cannot attain the CE 


in physically plausible setups I9l-ll2l. Another famous 
model is an information engine 13 which performs 
autonomous control and always transports a single par¬ 
ticle between particle baths. In contrast to Eeynman’s 
ratchet and the Biittiker-Landauer system, the informa¬ 
tion engine attains the CE. In addition, limiting effective 
filters of energy or chemical potential |lj-l2l| , a quantum 
dot [ 2 ^ [ 2 ^ , and soft nanomachines [ 2 J also attain the 
CE. However, contrary to externally-controlled engines, 
the comprehensive understanding of autonomous engines 
with finite difference of temperatures or chemical poten¬ 
tials has been elusive. Especially, the understanding of 



FIG. 1. (Color online): (a) Schematic of the autonomous 
Carnot engine, which consists of V-wall and H-wall. Since V- 
wall is at Z, particles can be exchanged with the particle bath 
with PH- (b) The potential landscape of V-wall. If H-wall 
is at t or b, V-wall is movable. Otherwise, V-wall is trapped 
at Z or r by delta-function type potentials, (c) Schematic of 
the torque. The red disk is on the cross point of V-wall and 
H-wall, which slides on both walls. The yellow disk is fixed 
and serves as a shaft of the gray slat. With a single rotation 
A, the gray slat makes one rotation and 
we extract mechanical work. 


macroscopic autonomous engines is missing. 

In this Rapid Communication, we address the issue 
of the general condition for autonomous engines to at¬ 
tain the CE with finite chemical potential difference. To 
demonstrate this, we first introduce a schematic model; 
an autonomous version of the macroscopic Carnot engine. 
This model clearly illustrates the characteristics of au¬ 
tonomous engines that in normal nonsingular setups they 
cannot attain the CE even with infinitely slow dynamics. 
Contrary to this, in the case of singular transition rates, 
since this singularity prohibits the particle leakage from 
the dense bath to the dilute bath, this engine attains the 
CE. We then move to general discussion and prove that 
without a special type of singularity any autonomous en¬ 
gine cannot attain the CE, which is consistent with both 
our model and the existing models 

Kinetic model and its coarse-graining. —It is hard for 
autonomous engines to attain the CE, and they attain the 
CE only if they possess a special type of singularity. To 
demonstrate the above characteristics of autonomous en- 
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gines, we introduce an autonomous Carnot engine, which 
extracts mechanical work from particle flux between par¬ 
ticle baths with given chemical potentials and fiL. 
The engine consists of two movable walls; V-wall and H- 
wall (see Fig.[T](a)). Only when H-wall (V-wall) is at the 
position t or b {I or r), V-wall (H-wall) can move along 
the X {y) axis, and otherwise V-wall (H-wall) is fixed at 
I or r {t or b) (see Fig. [TJ(b)). Thus, the engine has four 
stable positions, {l,b), and (r, 6), which we 

denote by A, B, C, and D, respectively. When V-wall is 
at the position I (r), the engine can exchange particles 
with the bath with yn (^i)- Otherwise, the engine can¬ 
not exchange particles. The engine is under isothermal 
condition, and the dynamics of the walls and particles are 
stochastic (the explicit time-evolution equation is shown 
in the Supplemental Material). Since the cross point of 
two walls passes through the rectangular-shaped trajec¬ 
tory A, hy imposing torque on 

the cross point we extract mechanical work automatically 
(see Fig. 12(c) and Fig. IS.II) . 

We here adopt a well-used approximation that with 
given position of the walls and particle number the parti¬ 
cles are in equilibrium [2^ [2^. If the equilibration of par¬ 
ticles is much faster than the dynamics of the walls and 
the exchange of particles with the baths, the above ap¬ 
proximation is justified. This time separation leads to the 
coarse-grained description; the Markov jump processes 
with discrete states (V, n), where X G {A, B, C, D} rep¬ 
resents the position of the walls and n represents the 
particle number. We call this model as coarse-grained 
autonomous Carnot engine (CGACE). 

The energy difference from Ato B including the exter¬ 
nal force is denoted by Eab, and EBC^ Ecd, Ejja are 
defined in a similar manner. Then the work against the 
imposed force per one rotation A^B^C^D^A is 
written as Eab + Ebc + Ecd + Ed a ='■ bFtot > 0. Since 
the engine is under isothermal condition, the transition 
rates of V-wall (i.e., A D and B gg C) should satisfy 
the local detailed balance condition: 

= -P{Exx~+F{Vx-,n)-F{Vx,n)), 

B[X -G- X;n) 

( 1 ) 

where we defined A := D,B := C, C := B,D := 
A, and Exx- '■= —Ex-Xy and /3 = l/feeT is the 
inverse of the Boltzmann constant times temperature. 
P{X -G X~\n), Vxy and F{V,n) represent the transi¬ 
tion rate from {X,n) to (V“,n), the volume at X, and 
the Helmholtz free energy with volume V and particle 
number n, respectively. Note that the details of the tran¬ 
sition rates of the H-wall as P{{A,n) -G {B,n')) are not 
important in the following discussion. 

Maximum efficiency of CGACE. —To confirm the dif¬ 
ficulty for autonomous engines to attain the CE, we here 
derive the maximum efficiency of the CGACE with fixed 
fiH and hl- In the following, we investigate the condi¬ 
tion for the maximum efficiency, and then calculate the 
efficiency under this condition. First, the engine with the 
maximum efficiency should prevent two kinds of leakage. 



FIG. 2. (a) Graphs of GA(n) and Gsin). The area of “abed” 
(colored by gray) corresponds to the rhs of Eq. (O, which is 
the upper bound for Htot- The area of “pbqd” (surrounded 
by bold lines) corresponds to the lower bound for G^. (b) A 
graph of GA{n) with singular transition rates ([9]) and (IIOII . 
Ga shows almost discontinuous behavior, which allows 
and nf, as n\ — nf, = 0(Vb^^®). 


One is leakage of particles: If the dynamics of H-wall 
is much slow, V-wall moves between B gg C or A gg D 
many times and particles leak from the dense bath to the 
dilute bath without extracting work. The other is the 
leakage of energy: If the exchange of particles between 
the baths and the engine is much slow, the walls rotate 
obeying the external force as A^D^C^B^A 
and the work is lost. Hence, it is plausible that the max¬ 
imum efficiency is realized when the dynamics of V-wall 
is much slow than that of H-wall and particles, and we 
treat this situation in the following. 

The stationary distribution Pst{B,n), for example, is 
then calculated as 


Pst{B,n) 


Pst{A,B) 


g—P{F{VB,n)GEAB—UHn) 

Zab 


( 2 ) 


where Pst{A,B) represents the stationary probability 
at A or B, and Zab := -b 

^-p(F{VB,n)GEAB-UHn) jg normalization constant. We 
denote the stationary probability flux of A —>■ X~ with 
n particles by jx^x-{n) := Pst{X,n)P{X -G A“;n). 
Owing to the law of large numbers, the realized particle 
number when the transition X ^ X~ occurs is around 

n*x ■= argmax„ jx^jf-(«)- 

We now calculate the efficiency y := Wtot/C^, where 
Cfi represents the average consumption of chemical po¬ 
tential per a single rotation A^B^C^D^ A. The 
condition that the direction of the dominant dynamics of 
the walls is A^B^C^D^A leads to 


jA^D{n*A) <jD^A{n*jj), ( 3 ) 

jc^Bijic) <jB^c(n*D)- ( 4 ) 


Summing the logarithms of ([3]) and and using the 
local detailed balance condition we arrive at a key 
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inequality: 


Wtot < 


{^Hng - - 



GB{n)dn 


{^Hn*A - fJ-Lrio - 



GA{n)dn I . 


(5) 



Here, Gx{n) {X = A,B) is defined as 
Gx{n) := ^ (^F{Vx,n) - ^lnP{X ^ X-;n)^ , (6) 


and n*x satisfies GAi'n’^) = Gsirig) = and 

Gsin’^) = GA{n*jj) = Note that the second law of 
thermodynamics implies monotonic increase of Gx(n). 
The right-hand side (rhs) of (O corresponds to the area 
of “abed” (colored by gray) in Fig. [2j(a). In addition, 
Cfj, is evaluated as > (fin — ^^L)[n*g — n'^), whose rhs 
corresponds to the area of “pbqd” (surrounded by bold 
lines) in Fig. [21(a). If a finite constant a < -foo satisfies 


OGa 

dn 


< a, 


n 


OGb 

dn 


< a 


for any n, the inequality ([S]) implies 


(7) 


IFtot ^ « (1 - 
Gfj, ~ pAfi 


^ 1 — ^Carnot 5 


( 8 ) 


where we used GA{n) < aln(n/n^) -|- for n > n*jj 
and GB{n) > aln(n/n|j) -f for n < n*g. The in¬ 
equality ([SI) indicates that the maximum efficiency of the 
CGACE is strictly less than the CE. Especially, if the 
transition rates of V-wall obey the symmetric rule [s^ 
or the Arrhenius rule [3l| , the condition d?]) is equivalent 
to the condition for the thermodynamic function of the 
gas: sup^ y n • jdr? < a. We note that the ideal gas 
satisfies n ■ d'^Fjdv? = /? for any n and V. 

In the foregoing discussion, it was shown that the 
CGACE cannot attain the CE with normal transition 
rates. However, the CGACE attains the CE with the 
transition rate with a special type of singularity. We 
again assume that the dynamics of V-wall is much slow 
comparing to the H-wall and the particle exchange. We 
set the transition rates between A and D as 

. g—fl{F{VD,n) — EDA) 

P{A D]n)= (,-l3{F{VD,n)-EDA) ’ 

. ^-PF{VA,n) 

P{D A;n)= ^-/3(F{VD,n)-EDA) ’ 


with a constant k. The transition rates P{B —>■ C;n) 
and P{C B;n) are also written in a similar manner 
with the same k. We note that such transition rates are 
physically realizable (see the Supplemental Material). 

The crucial point of the form of P{A —>■ D; n) is that 
Ga shows an almost discontinuous jump from /i^ to fiH 


FIG. 3. (Color online): (a) An instance of the state space of 
an engine. Xi^h , X 2 ,h, Xs^h touches the particle bath with 
fiH, and Xi^L, X 2 ,L touches the particle bath with fiL- Ar¬ 
rows represent possible transitions, (b) Schematic of the (dis¬ 
cretized) Biittiker-Landauer system and its state space. Xi^H 
(Xi^L) represents the position of the Brownian particle with a 
hot (cold) bath. The whole state of the system is determined 
by the position X, the energy E, and the direction of the 
motion c G {-b, —}. 


(see Fig. 121(b)). This discontinuity leads to the diver¬ 
gence of ndGA/dn, and thus the left side of the in¬ 
equality Eq. dH) does not prohibit the attainability of the 
CE. We then properly set Va,Vb,Vc,Vd = 0(Vb) and 
Eab,Ebc,Ecd,Eda = 0{Vo) as satisfying n\- n\) = 
0(Vo^^^) and n*g — = 0(Vb^^^), which are negligible 

in thermodynamic limit. Under this setup, the efficiency 
is evaluated as 


7 ] := 


lUtot 

Cm 


> 1-0 



( 11 ) 


which indicates the attainability of the CE with thermo¬ 
dynamic limit Vb —b oo (detailed setups and calculations 
are discussed in the Supplemental Material). 

Necessary condition to attain the CE. —We now leave 
the specific model and go to the argument on general 
autonomous engines with finite chemical potential dif¬ 
ference, which include models in Refs. [sl-l^. We note 
that if the engine has continuous variables, we take their 
proper discretization. We also note that our argument for 
particle baths is easily extended to the case of thermal 
baths. By adding a new stochastic parameter if neces¬ 
sary (see the Supplemental Material), an autonomous en¬ 
gine satisfies the following two conditions: (a) The engine 
touches at most one particle bath at one moment, (b) 
The energy difference of the engine between two states (as 
Exx- in the CGACE) is independent of its particle num¬ 
ber n. We then divide the possible states of the engine 
into {Xi,H,"- -.Xn.H} and {Xi^le" ,Xm,L}, where a 
state Xi^H{Xi^i) is attached to the bath with finitJ-L) 
(see Fig. 121(a) or Fig. IS.II) . The state of the whole sys¬ 
tem is written as (A, n). In the case of the CGACE, the 
engine takes four possible states {A,B} = {Xx^HiX^^h} 
and {G,D} = {Ai_l,A 2 ,l}- In the case of the (dis¬ 
cretized) Biittiker-Landauer system [^-[^, [T^, Xh{Xb) 
corresponds to the position of the Brownian particle at¬ 
tached to a hot (cold) bath (see Fig. 121(b)). 

We here use the fact that a thermodynamic engine 
which attains the CE satisfies the detailed balance con- 
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dition [s^ l- which is a consequence of the widely be¬ 
lieved conjecture that thermodynamic engines with the 
CE move quasistaticall y a nd has recently been proved 
for Markovian systems |33j . Then the stationary distri¬ 
bution for n with given X reduces to the grand canoni¬ 
cal distribution. Hence, the average of particle number 
transported from Xi^u to Xj^^ per unit time is calculated 
as 




E 


n ■ 


( 12 ) 


where we defined /(n) := P{Xi^H 

Xj^L]n). By assuming the local detailed balance condi¬ 
tion for X, {n)x- L^Xi h is calculated in a similar man¬ 
ner: 


{n)xjx^Xi^H 


E 


n ■ f{n)e^^^^ 


(13) 


Here, we used the fact that the local detailed balance 
condition is written in a similar manner to Eq. (El), which 
follows from the condition (b). 

Let Vb be a typical system size. Since the amount of 
extracted work is of order Vb, the necessary condition for 
absence of particle leakage between Xi^n and Xj^l is 


TT {{'n)xi,H^Xj,L - {n)xj,L^Xi,H) = 0- (14) 

Vo 

However, by compareing Eqs. E2D and m , monotonic 
increase of in terms of n yields 

Tj-{{n)xi,H^Xj,L - {n)xj,L^Xi,H) > 0 , ( 15 ) 

and the equality holds only when /(Vbp)e^^^'4)P has a 
delta-function type singularity in terms of p := n/Vb such 
that: 


jx^.H^x,,, (Y^P) cx cx 5{p - p*l (16) 

Yop) cx cx S{p - p*)- (17) 

Here jx^x'{n) represents probability flux oi X X' 
with particle number n. Note that p is not the parti¬ 
cle density. The conditions ES and E3 have a clear 
physical meaning: Particle numbers with both transi¬ 
tions Xi^H —>■ Xj,L and Xj^L are always the 

same unique value Vbp* within o(Vb), which is the only 
way to prevent the leakage of particles. Since pn and 
Pl are fixed, this singularity appears in two ways: (i) 
P{Xi^H —>■ Xj^L;n) = 0 for all n G N except n = Vop*. 
(ii) 1/Vb • d/dplnf{Vop) shows such discontinuity that 


p'—^p* — OVo—^oo Vq Op 

lim lim -^Eln/(Vbp) 

p'—>p*+o Vb—Vb dp 


< PL, 

P = P' 

(18) 

p=p' 

(19) 


To attain the CE, all possible transitions Xi^n ^ 
need to satisfy (i) or (ii). Without such singularity, the 
engine cannot attain the CE as seen in Eq. ([8]). This is 
our main result. 

Various existing autonomous engines which attains the 
CE [ 1 , [13I - It^ adopt the way (i), where always a single 
particle or a particular value of energy is transported. In 
this Rapid Communication, we construct an autonomous 
engine with the CE which adopts the way (ii). In con¬ 
trast, present autonomous engines which cannot attain 
the CE [^El satisfy neither (i) nor (ii). In the case of 
the Biittiker-Landauer system, for example, the momen¬ 
tum distributions of the hot bath and the cold bath at 
the contact point of these two baths are different. Due 
to this difference, round-trips of the Brownian particle 
between two baths cause the energy transport from the 
hot bath to the cold bath in the form of kinetic energy, 
which implies finite dissipation. 

Concluding remarks. —In this Rapid Communication, 
we derived the necessary condition for autonomous en¬ 
gines to attain the CE. The key property is a special type 
of singularity as Eqs. EH) and E3, which implies that al¬ 
ways the same and unique amount of particle number or 
energy is transported between two states with different 
baths. Without such singularities, an autonomous en¬ 
gine never attains the CE. This result is consistent with 
existiM results on the specific models of autonomous en¬ 
gines [a-IrTj. 

Such singularities are realized only by the way of trans¬ 
ports of a single particle or a particular energy (i) or a 
special type of thermodynamic limit (ii). Previous mod¬ 
els with the CE adopt the way (i), and no concrete model 
with the CE has adopted the way (ii). In this Rapid 
Communication, we constructed a concrete example of 
the way (ii), and thus it was shown that both ways (i) 
and (ii) indeed lead to the CE. Here, we should emphasize 
that the way (ii) is different from the following intuitive 
picture: For the case of the macroscopic heat engine, 
by decreasing the external force and slowing down the 
speed of the cyclic process, the engine will reach the CE. 
In fact, this intuition is wrong for almost all macroscopic 
engines. Although thermal fluctuation is usually negligi¬ 
ble in macroscopic engines, in the case of infinitely slow 
speed thermal fluctuation should be taken into account 
even in macroscopic engines and the above intuition over¬ 
looks this inevitable fluctuation. This fluctuation causes 
the leakage of particles or energy between two baths. As 
evidence of this, the inequality ([8]) still holds for the case 
of the infinitely slow speed. 

Our results may give a new perspective on the physics 
of molecular motors. Molecular motors are also au¬ 
tonomous engines, and in most cases they work with fi¬ 
nite chemical potential difference. Our results impose 
physical restrictions on molecular motors: To attain the 
CE chemical heat engines should adopt the way of sin¬ 
gle particle transports (i) or thermodynamic limit (ii). 
These two options look similar to the two types of molec¬ 
ular motors; working solely (like a kinesin) or collectively 
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(like myosins in a muscle) [s^. It will be interesting if 
these two characteristics are connected. 
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Supplemental Material 

Time-evolution equation and dynamics of the walls 

We here write down the time-evolution equation of the walls explicitly. Since V-wall and H-wall have the same form 
of time-evolution equations, we treat only the dynamics of V-wall. The x-coordinate of the positions of V-wall, /, r, 
and *-th particle in the engine are denoted by V, xi, x^, and Xi, respectively. The y-coordinate of the positions of H- 
wall, t, b are denoted by Y, yt, yb, respectively. Then the stochastic dynamics of V-wall is written as the overdamped 
Langevin equation: 


dX ^ dUviX) dUi{X,Y) ^ac/ 2 (x„v) , [2^^ 

= Fe.t,v(V, V) - - — -^ + V y 

i 


(S.l) 


Here, 7 represents the viscous coefficient, Fext,v(W, V) represents the imposed external force on V-wall at {X,Y), 
Uv{X) represents the potential energy of V-wall, Ui{X,Y) and U 2 {xi,X) represent the interaction energy between 
V-wall and H-wall, and between V-wall and i-th particle, and represents the white Gaussian noise. Ui{X,Y) is 
given by 


Ui{X, Y) :={S{X — xi — e) + S(X — xt + e) + S{X — Xr — e) + S{X — Xr + e)) 

• ixivt -e<Y <yt + e) + xivb -e<Y < yb + e)) 

+ ((5(V -yt-e) + S{Y -yt + e) + S{Y -yb-e) + S{Y - yb + e)) 

• ixi^i — £<X<xi+e)+ xixr — e < X < Xr + e)) (S.2) 

with sufficiently small e, where x(H) takes one if the condition P is true, and takes zero otherwise. 

We next discuss the autonomous attachment and detachment of particle baths. Only when V-wall is at r {1), the 
delta-function type energy barrier between the engine and the particle bath with is removed and particles 

are exchanged between the engine and the bath. This protocol is realized by using potentials as 

UH{xi,X) :=6(xi e S'//)(1 - x{xi - e < X < xi + e)), (S.3) 

Ul{xi,X) :=5{xi e 5 'l)(1 - xixr - e < X < Xr + e)), (S.4) 

where Sh (Sl) is a border between the engine and the bath with yn (ml) and S{xi € Sh) takes infinity if Xi is on Sh 

and takes zero otherwise. Note that such a protocol is realizable without any external control parameter. By tuning 





Particle number is conserved 



Particle number is conserved 




FIG. S.l. Schematic of the cyclic process of the autonomous Carnot engine. Sh {Sl) is a gateway between the engine and 
the bath with fj,H {yx), which is open at A and B {C and D) in this schematic. Since Sh is boundary, the dimension of Sh is 
one dimension lower than that of the engine. 
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FIG. S.2. (a) State space of the whole system of autonomous Maxwell’s demon. When the engine is at r (Z), particles are 

exchanged with a particle bath with {^l)- (b) Addition of a new stochastic variable e € { 61 , 62 }, which corresponds to a 
lever at right and left respectively. 


the external force and potentials for walls, the engine moves as A^B^C^D^A automatically and performs 
work against the force (see Fig. EH). 


Method to adding a new stochastic variable to satisfy conditions (a) and (b) 

Autonomous engines discussed in Refs. d S [H [Ml, 1M3 themselves do not satisfy the condition (a) or (b). 
By taking an information engine [ll | as an example, we here demonstrate that how the addition of a new variable 
makes these engines satisfying the conditions (a) and (b). 

The information engine is an isothermal engine with two particle baths with chemical potential and fSL- The 
engine has two possible states x G {r, Z}, and the engine possesses at most one particle (see Fig. IS.2l (aB. The engine 
exchanges particles with a dense (dilute) particle bath when the state of the engine is r (Z). The whole state of the 
system is determined by {x,n). We set the transition rates as satisfying 


P((r, 0 )^(r,l)) 

P((r,l)^(r, 0 )) 

(S. 5 ) 

, P((Z, 0 )^(Z, 1 )) 

P((Z, 1 )^(Z, 0 )) 

(S. 6 ) 

P((r, 0 )^(Z, 0 )) 

(S. 7 ) 

P((r,l)^(Z,l)) 

(S. 8 ) 


with Eq ^ El, which violates the condition (b). Here, we normalized /3 to 1. 

We now introduce an additional stochastic variable e € { 61 , 62 } to the engine (see Fig. IS.2l fb')'). The state of the 
modified engine is determined by {x,e), and the state of the whole modified system is determined by (x,e,n). We 
then set the transition rates as 


P{{x, 

n, 61) 

—>■ (x, n, 

62)) = P((x, 71,62) {x, 

n, 61)) = const 

(S. 9 ) 

P{{x,0,ei) 

-Z- {x, 

1,61)) = 

-- P((x ,0 

,62) -)■ (x, 1,62)) 

=P((x, 0 ) 

-)■ 


(S.IO) 

P{{x,l,ei) 

-)■ (x, 

0,61)) = 

= P((x,l, 

,62) -)■ (x,0,62)) 

=P((x,l) 

-)■ 

(x, 0 )) 

(S.ll) 





61) —>■ (Z, 1 ,61)) = 

= 2 P((r,l) 

-)■ 

G,l)), 

(S. 12 ) 




P{{IA: 

61) (g1,6i)) = 

= 2 P((Z, 1 ) 


(Gl)), 

(S. 13 ) 




P((r, 0 , 

62) —>■ (Z, 0, 62)) = 

= 2 P((r, 0 ) 


G,0)), 

(S. 14 ) 




P((Z, 0 , 

62) (^0,62)) = 

= 2 P((Z, 0 ) 


(gO)), 

(S. 15 ) 

F((r, 0 , 6 i) 

— j, (Z, 

0,61)) = 

= P((Z, 0 , 

61) (r,0,6i)) = 

=0, 



(S. 16 ) 

1—1 

— >■ (Z, 

1,62)) = 

-PiilA, 

62) (^1,62)) = 

=0. 



(S. 17 ) 


for X = l,r, and 
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Adding a new variable e 


FIG. S.3. (a) Schematic of the Feynman’s ratchet and pawl [^. The ratchet is driven by the cold bath and the pawl (i.e., 

green triangle) is driven by the hot bath, (b) The state space of (discretized) Feynman’s ratchet and pawl 0 and that of 
modified one. In the modified model, the additional parameter e distinguishes whether the system is attached to the hot bath 
or the cold bath. 


It is easy to check that the modified engine satisfies the conditions (a) and (b). Note that this addition procedure 
is different from the additional variable discussed in Ref. [l^ . By setting the transition rates for e much larger than 
that for X and n, the coarse-graining of the fast variable e leads to the original system written as (x,n). 

The addition procedures for other models are in a similar manner. As an example, we briefly discuss the case of a 
discretized version of Feynman’s ratchet [ll| . The state of the engine Xi is determined by the position of the pawl, 
which is discretized in two states, and the angle of the ratchet (see Fig. IS.3|) . The pawl is driven by the hot bath, and 
the ratchet is driven by the cold bath. At any state, both the position of the pawl and the angle of the ratchet can 
change, which violates the condition (a). We here add the new variable e G {e//, e^}, which divides the state of the 
engine Xi into a hot bath state Xi^h := [Xi^eu) and a cold bath state Xi^^ '■= {Xi,ei,)- Then transitions between 
Xi^H ^ ^j,H {Xi^L ^ ^ 3 ,l) are associated with only the hot (cold) bath, and thus the modified engine satisfies the 
conditions (a) and (b). 

We finally make a remark on this procedure. In the case of two heat baths we should introduce another heat bath 
which gives the stochastic nature of e. Since the energy with (A, ejj) and that with (A, e/,) are the same, the third 
bath does not exchange any energy with hot or cold heat baths. 


Setting of CGACE which attains the CE 

We now write down the detailed setting of the CGACE which attains the CE. The transition rates of the walls are 
given by 


P{A^ D;n) 
P{D A; n) 
P{B C- n) 


P{C ^ B-n) 


k . (>-P{F{VD,n)-EDA) 

(S.18) 

^-pF(VA,n) _|_ f,-l3iF{VD,n)-EDA) ’ 

. g-0F(VA,n) 

(S.19) 

^-pF(VA,n) _|_ f,-l3iF(VD,n)-EDA) ’ 

^ . g^-P(.F{Vc,n)+EBc) 

(S.20) 

g-pF{VB,n) f,-l3(F{Vc,n)+EBc) ’ 

h . p-pF{VB,n) 

f,-pF{VB,n) (,-l3{FiVc,n)+EBc) ' 

(S.21) 


Such transition rates are realized by setting the potential of V-wall as Eig. IS.41 The position A'{D') is close to A{D), 
where the distance of A — A' {D — D') are negligibly small. We set energy barriers with A — A' and D — D', 
large enough for V-wall to equilibrate between A' and D' before leaving to A or D. V-wall climbs the barriers with 
the Arrhenius rule: P{A —?► A'-,n) = 

We then set Va = Vb, Vb = riVo, Vc = rir 2 Vb, Vd = J'gVb, and take thermodynamic limit Vb —>■ oo with 
A 2 = 0 ( 1 ) and Eab,EbCi Ecd, Ejja = 0 ( Vb )- For sufficiently large Vb , the logarithms of the probability flux take 
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Equilibrate 

(Canonical distribution) 

FIG. S.4. The potential landscape of V-wall. Potential barriers with the height of exist at A — A' and D — D’. 

the following asymptotic forms with certain domains of n: 

InjA^oin) ~ - F{Va, n)) + ci 

^^3B^c{n) ~ j3{^iHn - F{VB,n) - Fab) + ci 
Injc-)-B(n) ~ l3{^iLn - F{Vc,n) + Ecd) + Cr 
\njD^A{n) ~ PiELn - F{VD,n)) + Cr 

Here, c/ := hi{kPst{l)/Z ab) and \Ti{kPst{r)/Z cd) are constants, and uad and ubc are the solutions of 

F{Va, uad) — F^VdtUad) + Eda =0, 

F{VB,nBc) — F(Vc,nBc) — Ebc =0, 

both of which are of 0{Vo). 

The relation (IS.22I) is confirmed as follows. The logarithm of jA^oin) is evaluated as 

InjA^nin) =ci - P{F{VA,n) - iinn) - ln(l + 

=ci - P{F(VA,n) - unn) - (S.28) 

Here, for sufficiently large Vq, f3{F{VD,n) — E{yA,n) — Ejja) with n close to but larger than uad is evaluated as 
P{F{VD,n) - F{VA,n) - Eda) = P{F{Vd , n) - F{VA,n) - F{VD,nAD) + F{Va, uad)) 


: n » hAD, 

(S.22) 

: n > UBC, 

(S.23) 

■ n < UBC, 

(S.24) 

■ n < UAD- 

(S.25) 

flBC solutions of 

(S.26) 


(S.27) 


~/3 


= 0{n - uad), 


dF 

dn 


{n - uad) 


1 


(S.29) 


whose sign is negative. Therefore, by setting n as n — fbAD = e^^^'FD,n)-F{VA,n)-EDA) i is satisfied, 

and since I3{F{Va, n) — fJ-Hn) = 0(Vb), the term is negligible in our analyses. Hence, we 

obtain Eq. (IS.221) . and in a similar manner we also obtain Eqs. (IS.231) . (IS.241) . and (IS.251) . We note that the fluctuation 
of the realized particle number around when the transition A ^ D occurs is = O{y/Vo), which is much 

smaller than 0(Vb^^^)- 

Now we write down all parameters in terms of Vq with fixed ri > 1 and (5 > 0 (the procedure is summarized in 
Fig. IS.51) . First, Eq. (IS.221) tells as the solution of 


dF 

dn 


= k'H- 


Vo,n\ 


To achieve n*D = n\ — under Eq. (IS.25|) . Vd is set as satisfying 


dF 

dn 


= k-L, 


VD,n\-Vo'^/^ 

where VdIVa ='■ ^2 = 0(1) holds. Owing to the extensiveness of the free energy, Eqs. (IS.23|) and (IS.241) yield 


(S.30) 


(S.31) 


dF 

dn 


VB,ri% 


dF 

dn 


Vc,n% 


= fJ-H 


= kL 


(S.32) 
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FIG. S.5. A flowchart of how the variables are determined by given Vo, ri, and <5. For example, using Vo = Vk, n\ is 
determined, and using Vo and n’X, Vd and n*o are determined. 


with rig = rin\ and n*Q = ring. Next, we set Eua and Esc as satisfying Eqs. (jS.26l) and (IS.271) with uad = 
{n\ + n*g)/2 and ubc = {n*B + ng)l2: 


Eda =F 



F 



nA + n*D 
2 


Ebc =F 



- F 



n-B + ric '^ 


(5.33) 

(5.34) 


This setting of Eda and Ebc ensures n\ > had > rig, n*g > ubc > rig, and \nx — rixx- \ = 0(Vo^^^), [rix- — 
nxx- \ = 0{Vo^^^), which confirm that n*x is close to the discontinuous region of G. Finally, we set Eab and Ego 
with a constant (5 > 0 as follows: 


Eab = F{yA,n*A) - F{VB,n*B) + ^JlH{n*B - n*A) - 5, (S.35) 

EgD = F(yc,ng)-F{yD,n*g) + ^lL{n*g-n*g)-5. (S.36) 


We now check that these settings satisfy Eqs. m and Owing to the stationary condition jD^A{n*g) — 

jA^D[ri*A) = jB^ c{n*B )-jc -^B{ri*g), it is enough to check both < jB^c{ri*B) anA jg^B{ri*g) < jD^A{n*g). 

Using Eqs. (IS.221) and (IS.231) . the logarithm of jB^c{ri*B)/jA^D{n*A) is calculated as 

^ _ F{VB,n*g) + F{VA,n\) - Eab) = <5 > 0. (S.37) 

jA^D[ri\) 

The logarithm of jD^Airig)/jc^Bin-g) is calculated in a similar manner. Hence, Eqs. ® and ® are satisfied. 
Finally, we evaluate Wtot and C^. First, VFtot is calculated as 


Wtot :=F{VA,n*A) - F{Vb, n^) + F{Vc, n*g) - F{VD,n*g) + HH{n*B - n*A) + ^lL{n*g - n*g) - 2S 
— F{Va, Nad) + F{VD,nAD) + F{VB,nBc) — F{Vc, ubc)- 


(S.38) 


Here, by using 


and 


we obtain 


OF 


dn 




OF 

dn 



n*A _ nAD _ 1 

^^0 ^0 " 2 Uo ^/ 3 ’ 


F{y^,n-A - = Vi (f(1, '&) - Fy^)) = v „. + ° (^)) = 


(S.39) 


(S.40) 


0(Uo^/^). (S.41) 
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In a similar manner, we obtain 


F{VB,nBc) - F{VB,n*B) = 

FiVc, n*c) - F{Vc, UBc) = 

F{VD,nAD) - F{VD,njj) = 0{Vo^/^). 

By combining n\ — n*j^ = and n*^ — tiq = and above four relations, Eq. (IS.38I) is evaluated as 

VEtot > W - [iL){n*B - n*o) + 0(Vo^/^). 


Next, Cfi is calculated as 


C,= 


ifJ-H - tJ-L)ie n*g + - n^) 

— 1 

= if^H - tJ-L)in*B - n*jj) + 0(Vb^^^)- 


0{e- 


Since n*g — n*jj = 0 {Vq), the efficiency is evaluated as 

VEtot 


V ■= 


Cu 


> 1-0 




(5.42) 

(5.43) 

(5.44) 


(S.45) 


(S.46) 


(S.47) 


which indicates attainability of the CE with Vq —>■ oo. 







